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ABSTRACT 

The precisions of extrasolar planet radius measurements are reaching the point at which meaningful 
and discriminatory comparisons with theoretical predictions can be made. However, care must be 
taken to account for selection effects in the transit surveys that detect the transiting planets for which 
radius measurements are possible. Here I identify one such selection effect, such that the number 
of planets with radius R p detected in a signal-to-noise ratio limited transit survey is oc R p , with 
a ~ 4 — 6. In the presence of a dispersion a in the intrinsic distribution of planet radii, this selection 
effect translates to bias b in the radii of observed planets. Detected planets are, on average, larger by 
a fractional amount b ~ a(cr / (R p )) 2 relative to the mean radius (Rp) of the underlying distribution. I 
argue that the intrinsic dispersion in planetary radii is likely to be in the range a = (0.05 — 0.13)i?j, 
where the lower bound is that expected theoretically solely from the variance in stellar insolation, and 
the upper bound is the 95% c.l. upper limit from the scatter in observed radii. Assuming an arbitrary 
but plausible value of <r/(R p ) ~ 10%, and thus b ~ 6%, I infer a mean intrinsic radius of close-in 
massive extrasolar planets of (R p ) = (1.03 ± 0.03)i?,/. This value reinforces the case for HD209458b 
having an anomalously large radius, and may be inconsistent with coreless models of irradiated giant 
planets. 

Subject headings: planetary systems - planets and satellites: general - stars:individual (HD 209458) 



1. INTRODUCTION 

After many years of effort, transit searches for 
extrasolar planets have recently come to fruition, 
with the detection of six extrasolar planets via 
the transit method. Most of these planets were 
discovered in deep Galactic field surveys by th e 
OGLE collaboration ljUdalski et al.l l2002alHEl [2003), 
and subsequently confirmed via radial velocity f ollow- 
up fiConackiet all l2003albt iKonacki et alJ l200l | 2005f 
Bouchv et all2004UPont et alJl2004UMoutou et al.l2004[ 



Pont et all 1200.1 [Bouchv ct al. 2005). One was dis- 
covered in a shallo w, wide-angle field survey by the 
TrES collaboration ijAlonso et all I2004D. Along with 
the transiting planet H D209458b l|Henrv et all 1200(1 
iCharbonneau et al"ll2000D. originally dis covered via ra- 
dial velocity survevs l|Mazeh et alJl2000|) . seven transit- 
ing planets are currently known, with many more sure to 
follow. 

The photometric light curve of a star with a transit- 
ing planet, when combined with detailed spectroscopic 
follow-up including precise radial velocities, yields the 
planet radius R p , mass M p , and semimajor axis a, along 
with the mass M*, radius and effective temperature 
of the host star. Of the routinely observable properties 
of transiting extrasolar planets, the radius is the primary 
diagnostic for testing theoretical models. Currently, the 
radii of the seven known transiting planets are deter- 
mined to precisions of AR p /R p = 5R P ~ (5 — 10)%. A 
number of different sources of uncertainty contribute to 
the total error, but ultimately the uncertainty in the ra- 
dius of a transiting planet is limited by the uncertainty 
in the primary mass, such that 8R P ~ (1/3)<5M*. Stellar 
masses can be inferred with precisions of ~ (10 — 20)%, 
and thus, in the future, we can expect to be able to mea- 

Electronic address: sgaudi@cfa.harvard.edu 



sure planetary radii to ~ (3 — 7)%. 

The steadily increasing number and precision of 
planet radius measurements, along with improvements 
in the theoretical models of irradiated giant plan- 
ets, has led many authors to present detailed com- 
parisons between measured radii and model predic- 
tions. Other than the mysterious case of HD209458b, 
whose anomalously large radius continues to elude 
a satisfying explanation llGui llot fc Showman! I2002t 
i Burrows et al.l l2003t iBaraffe et alJ l2003t iDeming et alJ 
20051 iWinn et al.ll2005|) . these comparisons have gener- 
ally res ulted in agreement betwe e n theory and obser- 
vations feodenheimer et all 120031 IBurrows et alj 120041: 



IChabrier et al.ll2004HGuillot]l2005tlLaughlin et al 12003) . 

Although detailed inferences about the ensemble prop- 
erties of giant close-in planets are not yet possible, the 
expected precisions should be sufficient to distinguish be- 
tween, e.g., models with and without a large solid core. 

Before celebrating the agreement between the theoret- 
ical predictions and observational constraints on extra- 
solar planetary radii, we should be certain that all rele- 
vant biases in the observati ons and/or theoretical predic- 
tions have been identified. IBurrows et alJ l)2003fl identi- 
fied one such bias, such that the planet radius as inferred 
from a transit light curve is ~ (3 — 10)% larger than 
the photosphcric radius predicted by models, due to the 
fact that rays from the primary passing perpendicular 
to the planet radius vector suffer a longer path length 
through the atmosphere. The purpose of this Letter is to 
point out that there exists an additional selection effect 
on the radii of planets detected in signal-to-noise ratio 
(S/N) limited field transit surveys. This selection effect 
is such that the number of detected planets is oc R p , with 
a ~ 4 — 6. In the presence of intrinsic scatter in the dis- 
tribution of planetary radii, this translates directly into 
a bias in the mean radius of detected planets relative to 
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the intrinsic population. This bias can affect compar- 
isons of theoretically-predicted radii with the ensemble 
distribution of observed planet radii. It will also affect 
the interpretations of individual systems if there exist 
unaccounted-for dependences of the radius on unobserv- 
able parameters of the system, such as the migration 
timescale or an indeterminate amount of internal heat- 
ing. Note that th e bias identified here operates in the 
same sense as the IBurrows et al.l lj2003|) 'transit radius' 
effect. Therefore, these two biases can easily combine 
to yield ~ 10% differences from the usual naive compar- 
isons. These biases are easy to understand. However, 
they must be acknowledged and considered when draw- 
ing conclusions about the agreement (or lack thereof) 
between theory and observations. 

2. SELECTION EFFECTS IN FIELD TRANSIT SURVEYS 

Field surveys for transiting planets are subject to 
a number of selection effects, which can lead to ob- 
served distributions of planetary parameters that are 
biased with respect to the underlying intrinsic distri- 
butions. These selection effects have been discussed 
bv lGaudi et~aT1 <|2005() and iPont et"afl l|2005j) . although 
these studies were primarily concerned with biases in the 
periods of detected extrasolar planets and did not con- 
sider the_bjaswrth respe ct to planetary radius in any de- 
tail. iPepper et all l)2003|) and lGaudi et all l)2005|) present 
simple scaling relations for the number of planets de- 
tected in a S /N- limited field transit survey, as a function 
of the parameters of the star and planet. I review the 
basic steps here, but refer the reader to these papers for 
a more detailed derivation of these relations. 

The total S/N of the photometric detection of a plan- 
etary transit is, 

I - ^ (£) < 

where N tr = (R^/TTa)N to t is the total number of points 
in transit, N tot is the total number of data points in the 
light curve, 5 = (R p /R*) 2 is the depth of the transit, 
and <7ph is the single-measurement relative photometric 
precision. In terms of fundamental parameters of the 
planet and star, the S/N scales as 

| oc R^ 2 T! qR y/ 2 d-\ (2) 

where T eq is the equilibrium temperature of the planet, 
and I have assumed that a p h is limited by Poisson noise, 

— 1/2 

such that tjph ocL, d, where and d are the lumi- 
nosity and distance of the host star, respectively. 

I note that the various factors in equation @ are un- 
likely to be independent. For example, it may be that 
planet properties are correlated with host star mass and 
thus i?» . In addition, models of irradiated planets predict 
that their radii will likely depend on the amount of stellar 
insolation absorbed by the planet. Thus R p is expected 
to be correlated with T oq , although precise predictions 
for this correlation are hampered by an incomplete un- 
derstanding of the physical processes involved with the 
absorption and redistribution of the energy from the in- 
cident stellar radiation. There is also observational evi- 
dence that the mass of close-in giant planets is correlated 
with a, such that planets closer to their parent star are 



more massive llOaudi et all2005HMazeh et al!2005lL All 

else equal, this would lead to a correlation of Rp with a. 
lLaughlin et alJ l)2005f) argue that such a correlation be- 
tween R p and a largely explains the apparent discrepancy 
between the period distributions inferred from radial ve- 
locity surveys and transit surveys. This is unlikely to be 
correct, first because the evidence for a correlation of R p 
with a is weak, and second be cause, even for the stre ngth 
of the correlation quoted bv lLaughlin et ahl l|2005|) . the 
magnitude of resulting selection effect would be less than 
half that due to the inevitabl e selection effect ar ising from 
the direct dependence on a lIGaudi et aljl200l . 

Although these correlations are interesting topics for 
future study, their contribution to the selection effects 
are unlikely to dominate over the direct R p term in equa- 
tion (J2J. I will therefore ignore the other terms for the 
remainder of the discussion, and assume S/N oc R^d -1 . 

At a limiting (S/N) m i n , the maximum stellar distance 
out to which a transiting planet can be detected is there- 
fore <i max oc R p . The number of planets detected is 
proportional to the total volume over which a planet 
gives rise to a transit with S/N > (S/N) m ; n . Assum- 
ing a constant volume density n of stars and no dust, 
the number of planets that can be detected is therefore 
oc dj^ ax oc R p . If the intrinsic distribution of planetary 
radii is fi(R p ) = dn/dR p , the observed distribution of 
planetary radii f {R p ) will be 

f (Rp) oc h{Rp)R^ (3) 

with a = 6. I leave a as a free parameter to allow for 
departures from the ideal discussed below. 

There are a number of assumptions that enter into the 
derivation of equation © , some of which I h ave explicitly 
stated , and others of which are reviewed in lGaudi et alJ 
(2005). The most relevant for the present discussion are 
the assumptions of a constant volume density of stars, 
no dust, Poisson-noise limited photometry, and a S/N- 
limited survey. All of these are violated to some extent in 
the actual field surveys that have detected transiting ex- 
trasolar planets. S. Dorsher, A. Gould, & B. S. Gaudi (in 
preparation) calculate the expected scaling of the num- 
ber of detected planets with R p for the OGLE surveys, 
using a realistic Galactic model of the source star and 
dust distribution, including the joint distribution of host 
star luminosities and radii, and considering the actual er- 
ror properties of the OGLE photometry. They find that 
the number of detected planets is a high power of R p , 
but with a somewhat smaller exponent than a = 6. The 
range is a ~ 4 — 6, with lower indices expected for plan- 
ets with smaller a. The properties of the TrES survey 
have not been described in detail, so I will assume for 
simplicity that they are similar to the OGLE survey. 

3. THE RADIUS BIAS 
It is clear from equation © that due to the R™ selec- 
tion effect in S/N- limited transit surveys, the observed 
distribution of planet radii will be a biased subset of the 
underlying intrinsic distribution. For example, the mean 
of the observed distribution will generally be larger than 
the mean of the underlying distribution. As a result, 
if this bias is not taken into account, incorrect infer- 
ences about the population of planets as a whole could 
be drawn from the properties of observed planets. 
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Fig. 1. — Bias vs. intrinsic dispersion in radii. The solid 
curve shows the fractional difference between the observed and 
intrinsic mean radius (the bias) versus the intrinsic dispersion in 
radii, assuming the intrinsic distribution is a Gaussian. The dotted 
curve shows a simple approximation derived from equation CI . 
The long-dashed curve shows the bias assuming there exists an 
additional population of planets whose radii are inflated by 10%, 
and constitute ~ 10% of the total population of planets. The 
shaded region shows the likely range of intrinsic dispersion in radii. 



The magnitude of the bias will clearly depend on the 
intrinsic distribution of planets. To provide a quantita- 
tive estimate of the bias, I adopt a Gaussian form for the 
intrinsic distribution /; of planet radii, with mean (Rp) 
and standard deviation a. The mean (R p } and stan- 
dard deviation a Q of the observed planet distribution f a 
can then be calculated in the usual way. Although exact 
analytic expressions for (R P )o and a can be found, they 
are cumbersome and not particularly illuminating. The 
bias, defined here as the fractional difference between the 
mean of the observed distribution relative to the mean 
of the intrinsic distribution, b = ((R p ) — (Rp))/(Rp)i is 
shown in Figure ^ F° r a /(Rp) ^ 1? the observed dis- 
tribution f is approximately a Gaussian. Therefore, the 
mean of f a can be approximated by its maximum, which 
yields, 



{R p )o — (Rp) 



1 + i ( 1 

2 2 



4a- 



1/2' 



(4) 



The resulting bias using this approximation is shown 
in Figure ^ it agrees well with the exact result. For 
a/(R p ) « (4a)- 1 / 2 , 



( a 



\(R P ) 



(5) 



Thus, if the intrinsic distribution of planet radii has a 
standard deviation of 10%, the mean of the observed 
planets will be larger than that of the underlying pop- 
ulation by ~ 4 — 6%. I note that the RS selection ef- 
fect also tends to result in an observed dispersion that 
is smaller than the intrinsic dispersion. However, for 
a I (Rp) 25%, the intrinsic and observed dispersions 
differ by < 10% for a < 6. 



Fig. 2. — Radius versus mass for the seven known transiting ex- 
trasolar planets. The circles show planets confirmed from follow-up 
observations of candidates from the OGLE deep field survey, while 
the square shows TrES-1, which was detected from a wide-angle, 
shallow field survey. The triangle shows HD209458b, which was 
originally discovered using radial velocities, and as such is not sub- 
ject to the same biases as the planets detect ed in field surveys. 
The s haded regions show the predictions of Bodcnhcimcr ct al. 
(2003) for radii of irradiated planets with equilibrium temper- 
atures in the range expected for the known transiting planets, 
T cq = 900 — 1800 K. These radii have been augmented by 
5%, to roughly accou nt for the 'transit radius' effect discussed in 
IBurrows et al. (2003). The dashed line and hatched region shows 
the mean radius and error of the intrinsic population of extrasolar 
planets, as inferred from 'debiasing' the data from the six planets 
detected in field surveys, assuming an intrinsic dispersion in radii 
of 10%. 



The magnitude of the bias is relatively insensitive to 
the form of the intrinsic distribution fi(R p ). As an ex- 
ample, consider the case where bulk of planets have the 
form assumed in the above analysis, but there exists an 
additional population of "inflated" planets. The inflated 
planets have the same dispersion as the other planets, 
but their mean radii are larger by 10%, and are ten times 
less common. Such a population might be indicated by 
the existence of HD209458b, which appears to have an 
anomalously large radius. In this case, the bias is very 
similar, as shown in Figure ^ 

4. SUMMARY AND DISCUSSION 

I have argued that there exists a bias in the radii of 
extrasolar planets detected in field transit surveys. This 
bias arises from the fact that there exists a strong se- 
lection effect in these surveys, such that the number of 
planets detected is oc R^, with a ~ 4 — 6. The mag- 
nitude of this bias depends on the form of the intrinsic 
distribution of planet radii, and is b ~ a(a / (R p )) 2 for 
a Gaussian distribution with mean (Rp) and dispersion 
a. Although the assumption of a Gaussian distribution 
of radii is reasonable, I stress that a bias exists for any 
intrinsic distribution with a finite dispersion. 

One might be tempted to conclude that this bias only 
effects inferences about the ensemble distribution of plan- 
etary radii, and that point-to-point comparisons between 
calculations and measurements of the radii of individual 
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systems would be immune to bias. This is not necessarily 
the case, however, because there may exist unobservable 
or poorly constrained parameters that effect the plane- 
tary radius, such as the migration timescale, age of the 
system, or an indeterminate amount of internal heating. 
These hidden parameter dependences may give rise to a 
dispersion in planetary radii at fixed values of the ob- 
servable parameters. Thus an observed planet is likely 
to have a larger radius than might be expected based 
on models incorporating typical values of these hidden 
parameters. 

Although the primary purpose of this paper is to sim- 
ply point out the existence of a radius bias, it is in- 
teresting to reconsider the comparison between models 
and data in light of this previously unaccounted-for ef- 
fect. To do so, I use th e predictions of mo d els of irra- 
diated giant planets by iBodenheimer et all l)2003|) and 
lLaughlin et alJ l)2005f) . Predictions of other models are 
similar, although there are discrepancies at the ~ 10% 
level. I use these models simply because these authors 
present their predictions in an easily accessible numerical 
form, specifically as tables of predicted radii as a functio n 
of planet mass M p and T eq l|Bodenheimer et al.l [2 0031. 
and as predictions for specific systems ^Laughlin et al.1 
2005). These models span planet masses in the range 
M p = (0.11 — 3)Mj, equilibrium temperatures in the 
range T cq = (113 — 2000)K, and consist of models with 
and without solid cores of mass 20 M® for M p < Mj 
and 40 Mg for M p > Mj. I increase all model radii 
by 5% to roughly account for the 'transit radius' effect 
discussed bv iBurrows et all l)2003f) . This is somewhat 
larger than the magnitude expec ted for more massiv e 
planets such as OGLE-TR-56b ijBurrows et alJ 120041. 
and about half that expected for less m assive plan- 
ets like HD209458b llBurrows et al.l 120031) I a ugment 
the specific predictions of lLaughlin et all ((2005) to in- 
clude the recently-confirmed planet OG LE-TR-10, us- 
ing th e stellar and planet parameters from lKonacki et al.l 
(12005ft . but accounti ng for the improved photometry of 
Hol man et al.l (l200oTl. Linearly interpolating from Tabic 
1 of IBodenhe imer et al.l l|2003|) . the predictions for the 
radius of this planet (T cq = 1317 K) are R p = 1.04 Rj 
and R p = 1.18 Rj for models with and without a core, 
as compared to the (preliminary) measured radius of 
R p = (1.16±0.05)i?j. 

Figure |5] shows the radii of obser ved extrasolar plan- 
ets, to gether with the predictions of IBodenheimer et alJ 
(2003) for the range of T oq spanned by these systems, 
T oq ~ (900— 1800)K. As has been noted by numerous au- 
thors, the radius of HD209458b is anomalously large. It 
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is important to emphasize that this object is not subject 
to the bias discussed here, as it was originally discovered 
via radial velocity measurements. The average observed 
radius of the remaining six planets is 1.10-Rj, with an er- 
ror in the mean of 0.03i?j and an RMS of OMRj. The 
intrinsic dispersion is consistent with zero, however the 
true dispersion is likely considerabl y larger. Just from 
the va riance in T cq , the models of Bo denheimer et alJ 
( 2003) predict a dispersion of (3 — 5)%. Any additional 
variance in the properties of the systems that affect the 
planet radius (i.e. in the age, metallicity, or core mass) 
would give rise to a larger dispersion. The 95% c.l. up- 
per limit on the intrinsic dispersion is 0.13 Rj, respec- 
tively. Thus the intrinsic dispersion is likely in the range 
0.05 - 0.13 Rj. 

Without detailed knowledge of the intrinsic scatter in 
the radii of giant planets, it is difficult to assess the effect 
of the bias on the current dataset with any confidence. I 
will therefore proceed rather speculatively. For definitc- 
ness, I will adopt a = 0.1 Rj. Although arbitrary, this 
value seems plausible given expectations about the vari- 
ance in planet properties and histories. Assuming a = 6, 
this gives a bias of b ~ 6%. Using this to 'debias' the 
mean radius of the observed planets gives an estimate 
of the mean intrinsic radius of close-in massive extraso- 
lar planet s of (JQ = (1.03 ± 0.03)i ?j. I n contrast, the 
models of IBodenheimer et al.l l)2003l) and lLaughlin et all 
(|2005T) predict (R p ) = (1.16 ± 0.01)R. ; for coreless plan- 
ets with the observed masses and effective temperatures . 
Provided that the models of IBodenheimer et alJ ( 2003) 
span the full range of intrinsic planet properties, then 
this inferred value of the mean intrinsic radius implies 
that the majority of planets have a massive core. This 
also strengthens the case for HD209458b having a anoma- 
lously large radius. 

Although these results are tantalizing, I stress that the 
comparison between the measured radii and theoretical 
predictions presented here is only preliminary. A more 
detailed study should be performed using a full suite of 
theoretical predictions, and a more careful consideration 
of the observational and theoretical biases. 
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